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Abstract 

Some optimization problems coming from the Differential Geome- 
try, as for example, the minimal submanifolds problem and the har- 
monic maps problem are solved here via interior solutions of appro- 
priate multitime optimal control problems. 

Section 1 underlines some science domains where appear multitime 
optimal control problems. Section 2 (Section 3) recalls the multitime 
maximum principle for optimal control problems with multiple (curvi- 
linear) integral cost functionals and m-flow type constraint evolution. 
Section 4 shows that there exists a multitime maximum principle ap- 
proach of multitime variational calculus. Section 5 (Section 6) proves 
that the minimal submanifolds (harmonic maps) are optimal solutions 
of multitime evolution PDEs in an appropriate multitime optimal con- 
trol problem. Section 7 uses the multitime maximum principle to show 
that of all solids having a given surface area, the sphere is the one 
having the greatest volume. Section 8 studies the minimal area of a 
multitime linear flow as optimal control problem. Section 9 contains 
commentaries. 

Mathematics Subject Classification: 49K20, 58E20, 53C42. 
Key words: multitime maximum principle, minimal submani- 
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1 Origin of multitime optimal control 
problems 

The original results in this paper show that the minimal submanifolds 
(see [1], [3], [4], [22], [31]) and the harmonic maps (see [2]) are solutions 
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of multitime optimal control problems (see [5]- [30]). In this way we 
change the traditional geometrical viewpoint, looking at a minimal 
submanifold or at a harmonic map as solution in a multitime optimal 
control system via the multitime maximum principle. Problems of this 
type are among the most challenging in Differential Geometry and 
Control Theory while being among the most important for various 
applications. 

Near Differential Geometry, there are many domains of science con- 
taining multitime optimal control problems: Material Strength (the 
description of torsion of prismatic bars in the elastic case as well as 
in the elastic-plastic case), Fluid Mechanics (the motion of fluid sub- 
stances, Navier-Stokes PDE written as first order PDE), Magneto- 
hydrodynamics (Maxwell- Vlasov PDE, Navier-Stokes PDE) etc. Of 
course, to describe some m-dimensional objects as optimal evolution 
maps, a deeply understanding of the meaning of evolution is neces- 
sary. Our results confirm that the basic ideas of Lev S. Pontryaguin, 
Lawrence C. Evans and Jacques-Louis Lions are still surviving for 
optimal control problems governed by normal first order PDEs. 

2 Optimal Control Problem with 
Multiple Integral Cost Functional 

Let us analyze a multitime optimal control problem based on a multi- 
ple integral cost functional and m-flow type PDE constraints [14]- [15], 
[18], [19], [28], [30]: 

max /(«(•)) = / L(t,x(t),u(t))dt 
subject to 

ox 1 

Q^it) = ^aC*)^*). = -,n;a = l,...,m, 

u(t) eu, te ^ot ; x (°) = x o, x(t ) = x to . 

Mathematical data: t = (t a ) = (t 1 ,...,*" 1 ) <G i?+ is the multitime 
(multi-parameter of evolution); dt = dt 1 A ... A dt m is the volume 
element in RT:; £lot is the parallelepiped fixed by the diagonal op- 
posite points = (0, ...,0) and to = (tji ^o 1 ) w ^ich is equivalent 
to the closed interval < t < to via the product order on 
x : Qot — > R n , x(t) = (£*(£)) is a C 2 state vector; u : $7oi U C 
R k , u(t) = (u a (t)), a = 1, k is a C 1 control vector, the running cost 
L(t,x(t),u(t)) is a C 1 nonautonomous Lagrangian; X a (t, x(t),u(t)) = 
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(X^it, x(t),u(t))) are C 1 vector fields satisfying the complete integra- 
bility conditions (m-flow type problem), i.e., DpX a = D a X/3 {D a is 
the total derivative operator) or 



fdX dX p \du* 



ex? dx a 



dt a dtfi ' 



where [XqjX^] means the bracket of vector fields. The complete inte- 
grability hypothesis constrains the set of all admissible controls (sat- 
isfying the complete integrability conditions) 

U = {u:R™^U\ D p X a = D a Xp) 

and the admissible states. 

To formulate the multitime maximum principle we need the control 
Hamiltonian 

H(t, x(t),u(t),p(t)) = L(t, x(t),u(t)) + p?(t)X l a (t, s(t), «(*))■ 

Theorem 1.1 (multitime maximum principle; necessary 
conditions). Suppose that the previous optimal control problem, with 
X,X l a of class C 1 , has an interior solution u(t) G IA which determines 
the m-sheet of state variable x(t). Then there exists a C 1 costate 
p(t) = (pf(t)) defined over S7ot such that the relations 



dpf OH 



and 



dapp'j (t) n (t)\gn ot = ^ (orthogonality or tangency), 
Q x 'j QH 

Op{t) = —(t,x(t),u(t),p(t)), Vt G n oto , x(0) = x 



H u a(t,x(t),u(t), P (t)) = o, vt g n oto 



hold. 

Remark 1.2 If the optimal control u(t) G U is not an interior 
point, then instead of critical point condition we have 

H((t,x(t),u(t),p(t)) = maxH(t,x(t),u,p(t)). 

Theorem 1.2 Multitime Global Maximality. The foregoing 
problem has a global solution (x(-),u(-)) if and only if the multiple 
integral functional Jq Mq H(t,x(t),u(t),p(i))dt is incave with respect to 
u. 
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3 Optimal Control Problem with 
Curvilinear Integral Cost Functional 



A multitime optimal control problem whose cost functional is the sum 
between a path independent curvilinear integral (mechanical work, 
circulation) and a function of the final event, and whose evolution 



PDE is an m-flow, has the form [14], [15]- [17], [19], [23], [24], [30] 



This problem requires the following data: the multitime (multi- 
parameter of evolution) t = (t a ) G ; an arbitrary C 1 curve Fo to 
joining the diagonal opposite points = (0, 0) and to = (^o> •••>*o l ) 
in the parallelepiped fioto = [0, to] (multitime interval) in BJ? endowed 
with the product order, a C 2 state vector x : flot — > R n , x(t) = (x*(t)); 
a C 1 control vector u : f2ot ~* U C R k , u(t) = (u a (t)), a = 1, 
a running cost L a (t, x(t),u(t))dt a as a nonautonomous closed (com- 
pletely integrable) Lagrangian 1-form, i.e., it satisfies DpL a = D a Lp 
(D a is the total derivative operator) or 



the terminal cost functional g(x(to)); the C 1 vector fields X a = (X l a ) 
satisfying the complete integrability conditions (m-flow type problem), 
i.e., DpX a = D a Xp or 



where [X a ,Xa] means the bracket of vector fields. Some of the previ- 
ous hypothesis select the set of all admissible controls (satisfying the 
complete integrability conditions, eventually, a.e.) 

U={u:R™^U\ DfiL a = D a Lp, DpX a = D a X p , a.e.} . 

The set U does not contain always the constant controls, but it contain 
sure controls which are continuous at the right (in the sense of product 
order). 

The previous PDE evolution system is equivalent to the path- 
independent curvilinear integral equation 




subject to Tj^W = X* a (t, x(t),u(t)), i = 1, n, a = 
u(t) eU, te n otQ , x(0) = Xq, x(t ) = x to . 



m 
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where 704 is an arbitrary piecewise C 1 curve joining the opposite di- 
agonal points and t of the parallelepiped f2ot = [0, t] C f&ot = [0' *o]- 

In the multitime optimal control problems with path independent 
integrals, it is enough to use increasing curves. 

Definition 2.1 A piecewise C 1 curve 7ot : s a = s Q (r), r G 
[ T 0) r i]i s ( T o) = 0, s(ti) = to is called increasing if the tangent vec- 
tor s = (s a ) satisfies s° > 0, with = only at isolated points. 

If we use the control Hamiltonian 1-form 

H a (t,x(t),u(t),p(t)) = L a (t,x(t),u(t)) +p i (t)X t a (t,x(t),u(t)), 

we can formulate the multitime maximum principle 

Theorem 2.1 Multitime Maximum Principle. Suppose that 
the previous problem, with L a ,X' l a of class C 1 , has an interior solution 
u(t) G U which determines the m-sheet of state variable x(t). Then, 
there exists a C 1 vector costate p(t) = (pi(t)) defined over f2ot such 
that the following relations hold: 

|p[(t) = -^(t,x(t),u(t), P {t)), vt g n oto , pj(to) = o, 

Q x j QH 

Ofi(t) = -g^(t,x(t),u(t),p(t)), Vt G n oto , x(0) = x , 
-^(t,x(t),u{t),P(t)) = 0, Vt G O oto . 

Remark 2.1 If the optimal control u(t) £ W is not an interior 
point, then instead of critical point condition we have 

H a ((t,x(t),u(t),p(t)) = max H a (t,x(t),u,p(t)). 

u 

Theorem 2.2 Multitime Global Maximality. The foregoing 
problem has a global solution (£(■), «(•)) if and only if the curvilinear 



integral functional I H a (t,x(t),u(t),p(t))dt a is incave with respect 
Jr oto 

to u. 



4 Multitime maximum principle 
approach of variational calculus 

In fact we show that the multitime maximum principle motivates the 
multitime Euler- Lagrange or Hamilton PDEs. 
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4.1 Case of multiple integral action 

Suppose that the evolution system is reduced to a completely inte- 
grable system 

— (t) = <(t), s(0) = xq, t € n otQ C 
and the functional is a multiple integral 

/(«(■))= / X(t,x(t),u(t))dt, (I) 

where Qot is the parallelepiped fixed by the diagonal opposite points 
= (0, 0) and to = (*o> ■•■j*o 1 )' * ne running cost X(t, x(t),u(t))dt is 
a Lagrangian m-form. 

The associated basic control problem leads necessarily to the mul- 
titime maximum principle. Therefore, to solve it we need the control 
Hamiltonian 

H(t,x,po,p,u) = X(t,x,u) + pfu % a 
and the adjoint PDEs 

^) = -^,x(t),u(t)). (AD J) 

Suppose the simplified multitime maximum principle is applicable 
dH__dX_ 7 7 _ _dX_ t _ i 

Suppose the function X is dependent on x (a strong condition!). We 
eliminate pf using the adjoint PDE. It follows the multitime Euler- 
Lagrange PDEs 

°X-D a (°£\=0. 

dx l \ dx l c 

4.2 Case of path independent curvilinear in- 
tegral action 

First, suppose that the evolution system is reduced to a completely 
integrable system 

— it) = <(i), s(0) = xo, t € tt oto C (PDE) 
and the functional is a path independent curvilinear integral 

J(u(-)) = I L p (t,x(t),u(t))dt^ u = (<), (J) 
Jr oto 
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where roj is an arbitrary piecewise C 1 curve joining the points and 
to, the running cost oj = Lp(t, x(t), u(t))dt^ is a closed (completely 
integrable) Lagrangian 1-form. 

The associated basic control problem leads necessarily to the mul- 
titime maximum principle. Therefore, to solve it we need the control 
Hamiltonian 1-form 

Hp(t,x,po,p,u) = Lp(t,x,u) +piU l p 

and the adjoint PDEs 

w {t) = -^ {t ' x{t)Mt)) - {ADJ) 

Suppose the simplified multitime maximum principle is applicable, 
i.e., 

We accept that the functions Lp are dependent on x (a strong condi- 
tion!). Then (AD J) shows that 

Pi{t)=Pi(0) - I ^(s,x(s),u(s))ds 13 , 
J Tot ® x 

where Tot is an arbitrary piecewise C 1 curve joining the points 0, t € 
Slot - From the foregoing last three relations, it follows 

-^(t,x(t),«(t)) = dJ Pi (0) - 5} ^( s , x ( s )Ms))ds X . 



If X° are functions of class C 2 , then applying the divergence operator 
d_ 

dtn'- 



D 1 = we find the multitime Euler- Lagrange PDEs 



OX 1 OXly 

5 Minimal submanifolds as 
optimal evolutions 

The minimal submanifolds are characterized by zero mean curvature. 
These become an area of intense mathematical and scientific study 
over the past 15 years, specifically in the areas of molecular engineer- 
ing and materials sciences due to their anticipated nanotechnology 
applications. The most extensive meeting ever held on the subject, 
in its 250-year history, was organized in 2001 at Clay Mathematics 
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Institute. In spite of all these efforts, the traditional thinking about 
minimality was not change. 

Recently, I gave a new approach to the theory of minimal surfaces 
[22], [29] changing the traditional approach into a new one based on 
solutions in a two-time optimal control system via the multitime max- 
imum principle. More, our recent studies [5]- [30] show the efficiency 
of approaching some classical problems using techniques of multitime 
optimal control or multitime modeling. 

Let Qq t be an interval fixed by the diagonal opposite points 0, r <G 
i?™. On an n-dimensional Riemannian manifold (M,gij) we introduce 
the multitime controlled dynamics 

^ (*)=<(*), t = l,...,n;a = l,...,m; (PDE) 

t = (t\...,t m )en 0T , x l (0)=xl x 1 (t) = x\, 

where u = (u a ) = (itjj : Qq t — > R mn represents two open-loop C l 
control vectors, linearly independent, eventually fixed on the boundary 
dQo T . The complete integrability conditions of the (PDE) system 
restrict the set of controls to 

f du* du l R ] 

w ={« = W = W|a#(*)=ol( t )}- 

A C 2 solution of {PDE) system is a submanifold (m-sheet) a : x l = 
x^t 1 ,...^" 1 ). 

The metric g a p{t) = gij(x(t))u t a (t)u J j3 (t) determines the volume 



^det(g a p) dt 



/Qot 

of the m-sheet x(t), t £ Qqt and this defines the functional 



V(u(-)) = - f ^det(g^) dt. (V) 



/n . 

Multitime optimal control problem of minimal submani- 
folds: 

max V(u(-)) 
«(•) 

subject to 

dx l 

— (t) = <(t), i = l,...,n; a = l,...,m, 

u(t) EU, t G r^Qr; x(0) = Xq, x(t) = Xl- 
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To solve the multitime optimal control problem of minimal sub- 
manifolds, we apply the multitime maximum principle. In general 
notations, we have 

x = (x l ), u = (u a ) = «), p = {p a ) = (pf), i = 1, ...,n; a = 1, ...,m, 

X a (x(t),u(t)) = u a (t), L(x(t),u(t)) = -\J&et{g a p) 
and the control Hamiltonian is 

H(x,u,p)= pfX l a (x,u) +p L(x,u). 

Taking po = 1, the adjoint dynamics says 

dpf dH 



dt a " dx r 



(ADJ) 



On the other hand, we have to maximize H(x,u,p) with respect to 
the control u, hence 

dH _ dL 

This necessary condition (critical point) is also sufficient since the 
foregoing Hamiltonian H(x,u,p) is an incave function with respect to 
u. Indeed, the Lagrangian L(x, u) is an incave function with respect to 
u, because it has any critical point, and the Hamiltonian H is obtained 
from L adding a linear term. 

Having in mind that x l a = «' a , we eliminate pf using the adjoint 
PDE. It follows the multitime Euler-Lagrange PDEs 

dL n ( dL . 

D Q [ — )= 



dx' 1 V dx l a 

or 

^L- Da l J^L d 9(h \ = o. 
dx i a \dg^ dxi J 

Summarizing, we obtain 

Theorem 5.1. A C 2 solution of the previous optimal control prob- 
lem is a solution of the boundary problem 

^~ Da { 9 ^ t|) = °' X(0) = X °' X(t0) = Xt » 

i.e., it is a minimal submanifold. 

The familiar relativistic free particle and the dual string, which 
provides a dynamical model of hadrons, are described by the PDEs in 
Theorem 5.1. 
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6 Harmonic maps as 
optimal evolutions 



Traditionally, the Harmonic maps are solutions to a natural geometric 
variational problem motivated by some fundamental ideas from differ- 
ential geometry, in particular geodesies, minimal surfaces, and har- 
monic functions. Harmonic maps are also closely related to nonlinear 
partial differential equations, holomorphic maps in several complex 
variables, the theory of stochastic processes, the nonlinear field the- 
ory in theoretical physics, and the theory of liquid crystals in materials 
science. 

Our recent papers change the traditional geometrical viewpoint, 
looking at a harmonic map as solution in a multitime optimal control 
system via the multitime maximum principle. 

Let f^or be a interval fixed by the diagonal opposite points 0, r € 
R+. Let h a/ 3 be a Riemannian metric on i?™. On the Riemannian 
manifold (M, g$j ) we introduce the multitime controlled dynamics 

dx l 

— (t)=u l a (t), i = l,...,n;a = l,...,m; (PDE) 

t = (t 1 ,...,t m )€n 0r , x\0)=x i , x\t)=x\, 

where u = (u a ) = (u l a ) : Slor — > R mn represents two open-loop C l 
control vectors, linearly independent, eventually fixed on the boundary 
dflo T . The complete integrability conditions of the (PDE) system, 
restrict the set of controls to 

w = {« = (««) = («J.)|^f(*) = ^(t)}. 

A C 2 solution of (PDE) system is a map (m-sheet) a : x l = x 1 ^ 1 , t m ) 
We introduce the energy density ^h a/3 (t)gij(x(t))u l a (t)u :> l3 (t) of the 
m-sheet x(t), t € f^or and the energy functional (elastic deformation 
energy) 

E«)) = -\( h^(t) 9ij (x(t))<(tH(t) dt. (E) 

The energy density is defined by the trace of the induced metric g a p = 
gijU l a u^ with respect to the metric h a/3 . 

Multitime optimal control problem of harmonic maps: 

max E(u(-)) 
«(•) 

subject to 



10 



— (t) = u % a (t), i = l,...,n; a = l,...,m, 

u(t) 6W, t G ^or! ar(O) = Xo, x(r) = a?i. 

To solve the previous problem we apply the multitime maximum 
principle. In general notations, we have 

x = (x l ), u = (u a ) = «), p = (p a ) = (pf),i = l,...,n; a = l,...,m, 

X a (x(t),u(t)) = u a (t), L(x(t),u(t)) = --/i a ^(*) ftj (s(t))«j,(t)«J(i) 
and the control Hamiltonian is 

H(x,u,p) =PiX t a (x,u)+p L(x,u). 

Taking po = 1, the adjoint dynamics says 

dpf dH 



dt a ~ dx i 



(AD J) 



On the other hand, we have to maximize H(x,u,p) with respect to 
the control u, hence 

du l a du l a 1 

This necessary condition (critical point) is also sufficient since the 
Hamiltonian H(x,u,p) is a concave (and hence incave) function with 
respect to u. In fact, the Lagrangian L(x, u) is an incave function with 
respect to u, because it has any critical point, and the Hamiltonian H 
is obtained from L adding a linear term. 

We eliminate pf using the adjoint PDE. It follows the multitime 
Euler-Lagrange PDEs 

dL „ / dL\ 

D a tt- = 0. 



dx l \ dx^ 

Summarizing, we obtain 

Theorem 6.2. A C 2 solution of the previous optimal control prob- 
lem is a harmonic map. 

Remark. It is possible to extend the notion of harmonic maps to 
much less regular maps, which belong to the Sobolev space H l (N, M) 
of maps from (N, h) into (M, g) with finite energy. The above equa- 
tion is true but only in the distribution sense and we speak of weakly 
harmonic maps. 
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7 Minimal volume at constant area as 
optimal control problem 

Suppose that D is a compact set of R m = {(t 1 , t m )} with a piece- 
wise smooth (m— l)-dimensional boundary dD. The volume / dt l ...di 

Jd 

of the domain D is related to the flux of the position vector t = (t a ) 
through the closed hypersurface dD by the Gauss-Ostragradski for- 
mula 

ml dt\..dt m = f 8 a Rt a N p da, 
Jd JdD 

where TV = (N^) is the exterior unit normal vector field on dD. On the 

other hand, the area of dD is / da. Introducing a parametrization 

JdD 

on dD, whose domain is U C R m ~ l , we have da = ||jV||dr/, where 
TV = ||jV||TV and 77 is an (m — l)-form. 

Let us show that of all solids having a given surface area, the 
sphere is the one having the greatest volume. To prove this statement, 
we take the normal vector field M as a control - a very interesting 
idea of my PhD student Andreea Bejenaru - and we formulate the 
multitime optimal control problem with isoperimetric constraint 

max/ 8 a pt a NP{t)dn subject to j yj 1 8 ali M a {t)N^{t) drj = const. 

Using the Hamiltonian 

H = 5 al3 t a M P - p \]8 a pN a NP, p = const., 
the critical point condition, in the multitime maximum principle, gives 

Since the Hamiltonian is a concave function of J\f, the critical point 
is a maximum point. This confirms that D is the sphere ||t|| 2 < p 2 in 
R m . 

Remark. In the optimal control problems, the Stokes theorem 
reads: Let oj be a controlled (p — 1) - form, 1 < p < m, with compact 
support on the p - dimensional submanifold D. If dD denotes the 
boundary of D with its induced orientation, then J 9d uj = J D du>, 
where d is the exterior derivative. Consequently, we can use as an 
action functional either those defined by the left hand member of the 
Stokes formula or those defined by the right hand member, eventually 
with new constraints. 
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8 Maximal area constrained by 
a multitime linear flow 

Let (M, g) be a Riemannian manifold and (TM, G = g + g) be its 
tangent bundle. Let S : J2ot — > TM be an m-sheet and suppose that S 
is expressed locally by x % = x l (t),y l = y l {t) with respect to the induced 
coordinates (x l ,y' t ) in TM, and t £ Qot as the multitime evolution 
parameter. Then the m-sheet S = itS in M is called the projection 
of the m-sheet S and is represented locally by x l = x l (t). Now let 
£ = f be a vector field on M tangent to S. Let £ y = 

be the vertical lift to TM. Let y = be the Liouville vector field 

on TM. 

Let = £(x(i)), y(t)) be the solution of the linear controlled 
m-flow 

^(t) = A) a (t)e(t) + Ai +ja (t)yi(t) + Bj(tK(t) 

on TM. On the tangent bundle we use the area 1-form 

w = g (ftiWWV^) - <7 ij (x)y i (x)^'(x)) 

fy* = dy* + rj fc ^'dx fc , <5f = df + rj-fc^'da:*, 
where is the Riemannian connection on M. Giving a closed curve 
C on the image of (£(t),y(t)), we introduce the area 

a = \Sc {diMC^W - gijWWW) ■ 

To find this area, we introduce the pullback 

\ (9iMt))C(t)yi(t) - giMtWitHitj) dt a . 

Then, the curvilinear integral 

*mo) = \J c (9iMt))e(t)yi(t) - 9 iMwmm) 

is the area of a piece from the m-surface (£(t),y(t)) bounded by the 
curve C = n(C). 

We formulate the multitime optimal control problem: max M (.) <j(u(-)) 
subject to the foregoing controlled m-flow. The Hamiltonian 1-form 
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+Pi(t) (A) a (t)x>(t) + K +ja {t)y\t) + Bl(t)u a a {t)) 

+«(t) (A^\t)x\t) + AlX) a {t)yi{t) + B^\t)u a a {t)) , 

is linearly affine with respect to the control variables, i.e., this Hamil- 
tonian 1-form can be written as 

H a = L a + M a u a a , 

where 

M a (t)=p i (t)B i a {t)+q i (t)B^ +i (t) 

are the switching functions. In general, there will be no extremum 
unless control variables are bounded, in which case they are expected 
to be at the boundary of the admissible region (see, linear optimiza- 
tion, simplex method). Suppose — 1 < u a a < 1. When the multitime 
maximum principle is applied to this type of problems, the optimal 
control u*^ must satisfies 

1 for M a (t) > 
? for M a {t) = 
-1 for M a (t) < 0, 

for each a = 1, ...,m. This optimal control is discontinuous since each 
component jumps from a minimum to a maximum and vice versa in re- 
sponse to each change in the sign of each M a (t) (switching functions). 
The optimal control u* a is called a bang bang control. 

9 Conclusions 

This paper refers to basic problems in control the partial differential 
equations of differential geometry. The original results are meaningful 
and useful for explaining many real world phenomena based on opti- 
mal controlled multitime evolutions. They shows that some dreams 
issuing from the papers of Lev S. Pontryaguin, Lawrence C. Evans 
and Jacques-Louis Lions are now partially covered by the multitime 
maximum principle. Of course, to pass from the previous local the- 
ory to a global one, we need another more flexible formulation of a 
smooth multitime optimal control problem involving two Riemannian 
manifolds. 
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